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ABSTRACT 
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systes design is presented. The emphasis is on 
compensator design but the method is applicable to the 
design cf any system with free parameters. Signal 
representation and system response are in the time 


domain. 


Tke inputs reguired from the engineer are the 
systen ccnofiguration, the desired output response and 
the free farameters. A parameter Dinimization routine 
is ther used tc minimize a specific cost function and 
to set the free parameters. A graphical cutput of the 
desired response and actual system response is tken 
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The past two decades have been witness to an ever 
increasing use cf the digital computer. Engineering usage 
and design are prcetably the most important justificaticn for 
the large memcry, high speed computers of tcday. Classes in 
ccemputer applicaticn to engineering problems have teéccne 
part of tke established curricula at post universities. 
Numerous papers have appeared on the apflication of the 
computer to engineering problems [1] which until recently 
were unsclvackle or at best required long, tedious 


precedures wkich aqave only approximate soluticns. 


Contrcl system engineering has relied increasingly on 
ccmputer sSinulaticn of large scale systems to verify that 
design specifications have been met and tc make 
mcdificaticns to a system even before producing a protctype. 
There are Epregrams available to simulate virtually any 
electrical cireuit, or control syStem, elther in transfer 
furction form cr as a system of first order differential 
equations. Cthers draw the Bode, Nichols cr Nyquist flots 
cf cpen and clcsed locp systems. Some programs help design 
cc&gfensatcrs which use an iterative methcd to achieve the 
desired frequency response [2] [3]- Researchers continue to 


better adapt the computer to engineering usage. 


Cantalapiedra [4] has used an iterative method tc find 
the optizum reduced order sodel for large order systems. 
MacNamara [5] went even further and used an iterative method 
tc find the cptimum compensation for an aircraft autofilot. 
It wculd eprear that these technigues cculd ke extendec and 


applied to the direct simulation and design of ccntrol 
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system circuits in the time dcmain. 


The irtert of this thesis was to develcp a user oriented 
program which cculd simulate a wide variety of ccntrol 
Systems and determine the values of the gains, feles and 
zeres necessary to prceduce a desired response. A convenient 
means of data card input was to be provided to specify (1) 
the contrcl system which was to be optimized and (2) the 
desired reepense. A locally available function ninimization 
reutine (ECXELX) was to be used to optimize the sinuvlated 
system's cutrut. 


fo simulate the system which is tc be opftinized, 
transfer furctions which are commonly encountered were 
reauced tc first crder linear differential equations. These 
egtations were then programmed so. that the transfer furction 
bliccks cculd be connected in an arbitrary fashion by data 
card infut. Several common nonlinear transfer blocks were 
alse provided. The program will simulate the system with 
the kncwo parameters and then allow all unknown or 
adjustable farameters to he fixed by the ccgputer 


optimizaticn routine to achieve the desired response. 
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Il. PRCGRAM DEVELOPMENT AND IMELEMENTATION 


et oe ee ae oe eee Te Se Se ae 2 2 oe ee eee > a > ee ae ae cows 


A. GENEFAL 


Any frogram which is to be of maximum tkenifit tc the 
user must kave a simple means of data input and an cutput 
which is e€asy to interpret and apply te the problem at hand. 
The input data should have a physical significance that does 
nct lose its relevance through the programmirg of nunerous 
egtuaticns. The program should be a readily usable tocl and 
not a proktler in itself. The intent of this thesis is to 
"present suck a program, Computer Automated Lesign of Systems 
(CADS), which is readily used £OLr Sinulation with 
Oftimizaticn cf the output. Optimum in the sense that the 
output respeonse cf a given system configuration is as near 


the desired respcnse as possible. 


Most cecntrol system design starts with a prcposed 
transfer function block schematic to achieve a desirec time 
dcmain response. CADS has the ccmmon transfer functions 
built into the prcegram and the input data can come directly 
frcm the frofrosed system schematic. The transfer tlocks 
which are availakle for system simulaticn are presented in 
Tarle I. Tkese blecks should be adequate, either seperately 
or in varicus cascade combinations, to represent sost 


ccentrol systems. 
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CF SYSBOL EQUATICN SOLVED 


Bare -P Og + Gea 


Bo -P 65+ G (6;+26;) 
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TAELE I Pregram Transfer Functicn Blocks 
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To use the program the engineer must know the systen 
configuraticn in transfer block form and the desired second 
order cutput response. If other than a second order 
response is desired, this may be easily specified but 
reguires scre knowledge of how to input the desired 
response. CALS will take the transfer block system, ccnnect 
it, ccmpare the given system response tc the desired 
response and set any free parameters to achieve the closest 


match possible of the system output to the desired restcnse. 


Figure 1 is representative of the type of system the 
prcegram car simulate and optimize. In figure 1, the 
parameters of the numbered transfer function blecks are 
knewn or fixed by equipment limitations. MTransfer blocks X, 
Y, and Z ccntain variable parameters which must be selected 
by the prcgram to make the system cutput rerrcduce a desired 


output furcticn as closely as fossible. 


- = 


FIGURE 1 Typical System to be Crtimized 
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E. MAIN FRCCEAM 


The MAIN program was developed to contrel the selection 
of the sutfurcticns used in the optimization frocess and to 
ccgpute tke cesired response ci the systems which were to be 
cetimized. A second order ster response is calculated fron 
equation (I) as the desired response and stcred in the array 
XCATA. 

G 
Go = see 25useue “tt (I) 
The MAIN is essentially a bookkeeping routine which controls 
the Frcecraa execution. The subfuncticn operaticrs it 


ccentrols are defined in the following secticns. 


C. ELAN 


The ccmacn simtlation routines LISA, DSL, ECAP, CS¥*P and 
INTEG were irvesticated in an attempt to adapt them tc _ the 
ceneral frotlem specified above. These prcgrams required 
either the irput cf the system's state variable equaticns, 
Lapiace trarsform eéquations or the node-transfer function 
pair for each system simulated. These programs are useful 
for simulaticn if all of a system's parameters have been 
specified. Eawever, each simulation is problem specific and 
these teéechniguesS are not readily integrated with cther 
sukfuncticn pfregrars. A sSimulaticn subproaqram whickt was 
Capatle cf accepting data card input tc simuiate a wide 
Variety cf centrel systems and of working with other 


existing subrrogrags had to te developed. 


LS) 


Subfurcticn PLANT was developed to simulate the systems 
which were tc ke optimized. It provides the versatility 
necessary to simulate numercus system configurations and is 
capable cf werking with a wminimizaticn recutine to cptimize 
the variakle ,arameters. The program reads the transfer 
furction blcck connections from data cards. This data is 
then used to autcmatically set up the system equatiors in 
state variakle format. These state variatle equations are 
sclved as first crder, ordinary differential equations by 
the Runge-Kttta-Gill forth order method. The programming 
was done in Fortran IV and all calculations are in dcuble 
precision. The capability of connecting the transfer 
furction biccks in any configuration (feed forward, 
feedback,etc.) .1S possible with this program. Simuiation 
flexibility is provided by having each Elcck within the 


system capakle of accepting an external forcing functicn. 


1, Eleck Data Card 


Several possible ways of inputting the data 
necessary tc simulate arbitrary system configuraticns were 
investigated. A means of defining a system configuration 
directly frem a bicck diagram schematic was chosen as _ the 
mest preferakrle method. USing this approach, each data card 
was designed to be directly related to a transfer furction 
of the syste#. This resulted in having direct access to the 
transfer function farameters for optimizaticn. 


Ic simulate the system, a data card is prepared for 
€ach transfer bleck in the schematic diagram. The data card 
ccctains a field cf numbers which specify the block ntmker, 
tke type cf transfer function ccntained within the tErlock, 
tke input ncde and the output node to which the Elcck is 
cconected and the values of any parameters associated with 


tre transfer function. The general fcrmat of the data 
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cards used tc input the system configuraticn is shown telcw. 


{1 Poa {20 140 }60 (Column Numter) 
ELKCCLD=EFEVV G P Z 
Where: | 


CC = Pcsiticn number of the block 
D = Type of Elock (number) 

EE 
VV = Cutfut ncde number 


Input ncede number 


G, EF, and Z are parameters cf the transfer functicn. 


The rumker of transfer blocks in the system and the 
Gata cards associated with each of the blecxs are read upon 
initial entry into the simulaticn routine. The program then 
connects the transfer function blocks in the proper order by 
ccf@paring the input node number of a block to the cutput 
ncde nugter cf every other block. Whenever these two 
numbers are equal, the blocks are knowr to ke connected and 
a flag is set equal to 1.0. The flags are then used to 
identify the input drives to each block. 


The input quantity to each block is called _ the 
CRIVE. The program was designed to ailow for multiple 
inputs to the system teing simulated. This flexibility was 
actieved by mgaking the input to each block equal to tte sum 
of the outputs of all blocks connected to the input node 
rlus any external forcing function (DRVIN) feeding the 
transfer tlcck. The input DRIVE to a Elock is determined as 


shcwn in €guation (II). 
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DRIVE(i) = )> > THA(3j) *FLAG(35,i) + DRVIN(i) (II) 


may) 4«6is the output of block j. FLAG(j,1) is 1.0 if block 
j is connected to Elock 1 and zero if it is not ccnnected. 
DRVIN(1i) is any external forcing functicn specified ry the 
user whick drives tlock i. DRVINS must be inserted in 
sukroutine ELANT as Fortran IV statements. The standard 
prcegram has [CHVIN(1) = 1.0 specified as a unit step infut to 
the system. Proklem III-E in the section Investigation of 
Prcgram EFerfcrmance demonstrates how multiple, time varying 


inputs are tc ke inserted in the fprogran. 


{Ike transfer function blocks available for system 
Slimulaticn are shown in Table [. These blccks were selected 
because cf tkeir ccmmon usage in the modeling process. They 
were aisc fcurd to be adequate either separately cr in 
Various cascade combinations to represent most ccntrol 


systems. 


fre transfer function equations fcr each tyre of 
bicck were written in state variable format and stored in an 
array nagred THA. The program reads a number from a block's 
Gata card which specifies the type of transfer furction 
asscciated witk the block. The system eguations are then 
seclved seéguentially by selecting the type of eqvation 
assccitaed with bicck number one, solving fcr its output and 
then sequencing tc block number two, etc. The integral 
equations are solved by a modified RKL fcurth order nrethod 
in the surfunction rcutines. Three integration routines 
were necessary tc store the irtermediate results obtained 
for those eguations involving a double integraticn. The 
fceur sutfunctions, RKLDE2, RKLDE3, CCPLX and RKILDE4 are 
called Ly ELANT. EKLDE2 is used for the integraticn of a 
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mpe two blicck. RKLDE3 is used for integration cf a type 
three blcck and to store intermediate quantities for 
Seutruncticn CCPLX. CCPLX calls RKLDE3 and RKLDE4 for 
integraticn cf a type four block. 


Nc freovision has peen made for the input of initial 
ccnditicre cr the integrators. Therefore, the integrations 
must start at time egual zero. The user must specify the 
step size tc be used for integration (DT) and the _ fefroklen 
run time (IF). To conserve computer time [CT should be made 
as large as fpossible. DT = TF/1000 is suggested as 
appropriate in most cases. Equally important is tnat TF be 
as small as fessible. Use of the program has shecwn that 
letting IF te greater than the transient response time of 
the systen is rarely justified. For preliminary analysis fF 
was kept tc cnly slightly lcnger than the time of the first 
undershoct fcr second order dceminant systems. Because of 
the complexity of subroutine PLANT a flcw chart of PLANT has 


keen includec as Appendix B. 


{he parameters of the system which are to be 
oftimized ky the bininization routine EOXPLX must _ be 
specified in PLANI. The minimization routine returns the 
trial values cf the variable parameters to ELANT in an array 
lakeled 'C'*. Regular Fortran IV statements equate the 


Variable parameters to the values in this array. If a 


pele-zerc fair were to he optimized, the fcllowing 
statements wceuld be inserted in PLANT: P(i) = C(1) 
and Z(i) = C(2). The location fcr the preceding statements 


is clearly indicated in the program listing for PLANT. 
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D. DESIBED KESPONSE ({XDATA) 


The criteria against which the system response will be 
ccmpared will vary according to the applicaticn cf the 
system. Since most design work on control systems is done 
cn the tkasis of second order dominance, the progras was 
written tc provice a second order step response with 


adjustable 6, W and gain as the basic criteria against 
n 


which the sinulated system will be compared. The desired 6, 
‘. and gain are read in as input data. The step response is 
then ccmputed in the main program by subfunction RKLITEC and 
stored in the array called XDATA. Any cther time domain 
Lesponse may be specified by the user by remcving the second 
Order step response equaticns and replacing them with the 
eguaticns of the desired response. An example is fpercvided 
by problem IITI-E ir Investigaticn of Program Performance. 
If the pregram is EFeing used for simulation only and nc data 
curve is desired, setting the input variatle LEAP = 1 will 
cause the pregram to bypass the computation cf the standard 
secend crder data equation. 


E. COSI FUNCTION (PERFORMANCE INLEX, FI) 


The achieved system response is ccmpared witk the 
desired resjcnse and the difference is the errcer. The 
FrOgram sé€arckes for the parameter settings which will 
ainimize this error. The cost functicn may be specified by 
the user tc weight the system outputs as desired in 


sukfuncticn FE. The default cost function cf the program is 


TE 


tke integral e€rrcer sguared. Jj = (Err)? dt An exampl 
: € 


of a weighted cost function is presented in Section III 
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F. MINIMIZATION RCUTINE 


1. seneral 


The free parameters are optimized to reduce the cost 
function ty the complex methcd of M. J. Box.{6] Hox's 
ccnstrained crtimization method has been programmed at the 
Naval Fostgraduate School by R. BR. Hilleary as a sukrcutine 
called ECXELX. This subroutine will find the minimum cf an 
acrkitrary function (cost function) subject to arbitrary 
explicit constraints and for implicit constraints. Explicit 
ccnstraints are defined as upper and lower bcunds on the 
free parameters. Implicit constraints may be arbitrary 
functicns of the free parameters (e.g. PP, <9 178) . 


Iwo function subprograms are used to evaiuate the 
cbjective functicn and implicit constraints, FE and KE 
respectively. The method BOXPLX uses te search fcr the 
vaiues of tke free parameters which minimize the cost 


function is explained in the ccmputer program listing. 


2. Exgdicit/Implicit Ccnstraints and 


= —- —_ = 


ECXELX searches a feasibility region (n-dimensicnal 
Space, where n = number of free parameters) defined Ey upper 
and lower Ecunds on the free parameters for a minimus cost 
function. The sgaller the region defined by these 
koundaries, the amore rapidly the program will converge to 
the optingum parameter settings. Good engineering judcement 


will ke necessary to keep the feasibility region as smail as 
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possible. Tke boundaries of the search regicn are read fron 
Gata cards Ey the main program as the upper bound (XU) and 
the lower bcturd (XL) for each free parameter. 


Implicit constraints may be any arkitrary furcticn 
of the free parameter desired. If an implicit constraint 
Such as tne product of a poie-zero pair must be less than 
scme numker is to Fe evaluated, it must be supplied by the 
user to tke subfunction KE. No implicit constraints were 


used in tktis thesis. 


The starting values cf the free parameters (XS) are 
read in Ey tke MAIN frem data cards. A good choice of 
Starting values wili dramatically reduce the time required 
fer optinizaticn. A preliminary root locus or Bode plot 
methcd of estimating the best values of the free variabies 


shculd ke accomplished whenever possible. 
G. GRAFHICAL OUTEUT 


Two sukroutines were written to provide for the 
Qraphical output of the desired response and the best system 
response achieved ty the optimization frocess. The user may 
select, by data card input, either subroutine PPL& waich 
prcevides a high speed printer flot or subroutine PIC which 
prevides a calcomp graph. Every fifth integraticn point 
stered in the arrays xXDATA and THAOUT ius gflotted. 
Subroutines EFIT and FIC call the subrcutines PLOTP and LCRAW 
respectively. FLOIP and DRAW are standard flotting routines 
at the NES computer facility and are not a part ci the 
Sinulaticn frogram. Figure 2 diagrams the information flow 


and data input to the progran. 
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The example problems presented below were used to aid in 
the develcrment cf the optimization pregram. The order of 
difficulty cf the problems progresses frcm a Simple text 
bock single variable, single input system tc a nultivariakble 
Operaticnal servo drive system which has multiple inputs and 
discrete level feedback. An example of how a schematic 
diagram representaticn of a system 1S prepared for infut to 
the progran is presented pricr tc considering the example 


problems. 


A. CATA INECT, AN EXAMPLE 


The Ward-leonard drive system [7] shown schematically in 
Figure 3 (a) has two variable parameters. The gain cf the 
amplifier and the tachometer feedback are available to 
adjust tte system's response. To simulate the system a 
blcck diacraz representation of the system is drawn as shown 


in Figure 3 (c) using the transfer blocks frem Table Il. 


(a) 
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DRVING)+ 


FIGCKRE 3 Ward - Leonard Speed Ccntrecl Systen 
Using Feedback. (A) Schematic Ciagran 
(B) Elock Diagram (C) Bicck Diagram 
Using CADS Blocks. 


Where G(1) = K/r P(1) = WIfr 
G(2) = Km/yRa 
EUS) ee P(3) = 0 
G(4) = -Kn 
G(5) = -Kt 
DRVIN({(1) = Er 
DRVIN (3) = -T, 

The nodes cf the block diagram are then nugbered 
eeruentially 1,2,.-.<.«,0. The Elccks between the ncd€és are 
also nugbered seguentially M2 eee N as shcwnr nye 
Mragure 3S (c). Data cards (N) are then prepared for each 


bleck which specify the block number, type cf bicck, the 
input ncde, the output node, and the parameters ccntained 


meetin tbe bicck. In Figure 3 (c), for example, biock 1 is 


a type two transfer block ccnnected Ltetween nodes 1 and 2. 
The data carc input for this block would be 


BLKO 12=0102 K/, 1/, 


The program reeds the data card input and connects the 
biccks Ey setting a FLAG = 1. whenever the input node 
number to a Elock is the Same as the output nede number of 
any other tlcck. The input to a transfer block is then 
determined tc be the sum of the outputs of all tklccks 
ccnonected tc the input node plus any external fercing 
function driving the input node. The program has a unit 
step specified for DRVIN (1). If this is the only infut to 
the system, ne action is necessary on the part of the user. 
If cther than a unit step input to the system is desired or 
if there are other external forcing functions such as 
CRVIN (3), they must te specified and placed within the kody 
cf the sukfurction PLANT as Fortran IV statements. Fer the 
example shcwn in Figure 3 (c), a card with the equation 

DRVIN(3) = £(T ) 


would have tc be inserted preceding the drive equaticns. An 
example cf hcw multiple, time varying drives are specified 


is given in section III. E,. 


When cptinizing a system, scme of the input quantities 
Will be unkrcwn or variables. These variables must also be 
assigned witkin subroutine PLANT. To optimize the variables 
2 and a cf Figure 3 the follcwing two statements would be 


inserted in ELANT: 


G (1) 
G (5) 


Ci} 
C (2) 


where C(1) and cC(2) are the variables which will be 


optimized by subroutine BOXPLX. 


Za 


E. TACHCMETER FEEDBACK 


The first optimization protlem attempted was one which 
had an exact solution that can be found by algékraic 
methods. Arn instrument servo [8] with unity feedkack and 
fcerward transfer function 

1000 
S(S+10) (III) 


was to be ccupensated with tachometer feedback as shcwn in 


G (S) 


Figure YW. The cnliy specification fcr the system's 
performance cf this single variable, seccnd crder system was 
the simple requirement that the closed lccf roots have a 
peee= 0.7. 


$@eo 
S$ (S+ 10) 


K,S 


FIGURE 4 Tachometer Compensated Systen 


The system shewn in Figure 4&4 was redrawn as Figure 5 in 
crder Ee achieve the tachometer feedback. The 


characteristic equation of the system shown in figure 5 is 
S2 + (10 + 103 le + 103 = 0 (IV) 


Tke required i = 0.0343 may be calculated from egtvaticn 
eV 
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FIGURE 5 CADS Block Diagram of Tachcreter Systen 


CADS was frogramazed to optimize the system with the 


variable, Re specified to be between the limits 
0.C1 < A < 1.0. The desired regponse was specified tc be 
the standarc seccnd order step respense for 6 = Q.7, 
Wn = 100C. CADS determined the optimum value for Se to be 
a = 0.0343. Figure 6 shows the system's step respense and 


the desired response are virtually superimpesed. 


FOS, 
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FIGUFE 6 CADS Compensated Tachcmeter Feedback 


System Response 
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The type one system, in the preceding example, allowed 
derivative feedback from the forward path without regqtiring 
a block which was capable of differentiation. One may 
enccunter a type zero system or a system which cannct be 
ccenfigured to previde derivative feedback from the fcrward 
path. The pessibility of providing a fseudo derivative 
feedback using a type three bleck was investigated fcr these 


cases. 


A type three bicck with Z = 0 and G = P is shcwn in 
Ficure 7. If the pole is placed far out or the real axis, 
this blilcck afrroximates derivative feedback. 


FIGUKE 7 Type Three Block used for 


Pseudo Tachcmeter Feedback 


The effect cf using this pseudo tachometer feedrack on 
the ccmplex roots of the closed loop system is neglicikle. 
It dees acd @ real roct at P = -964. 


The system of Figure 5 was redrawn as shown in Ficure 8 


using the pseudc tachometer feedback. 


FIGURE 8 Pseudo Tachometer Feedback 
Tre optinizaticn prcegram calculated Se =——0203592 £66 5ehe 


atcve configuration. This gave a 46 = 0.715. 


The system response and the desired response are shown 
woe Ligure S. The two responses are again rearly 
superinpcsed. This method cf providing derivative feedkack 
has the disadvantage of adding an integraticn step tc _ the 
Eroblerm sciuticn with the concomitant increase in prctlen 


ecluticn time. 


THAOUT 
0.4 


FIGURE 
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Time (sec) 


Fseudo Tachometer Ccmpensated System's Response 
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C. CASCADE CCMEENSATION 


Having jfreven the feasibility of the prcegram optinizing 
a single variable system where an exact solution was 
availakle, tke next yroblem considered extending the froblen 
sccpe tc a third order plant with two variables. The 
unstable plant that was tc be compensated is shcwn in 
Figure 10. 


—————— 
Cjul(jwrC2jurn 


FIGUKFE 10 Third Order Plant to be Compensated 


The slant was te be stabilized using a single section 
cascade ccmpensatcr. The compensated plant was reguired to 


have - <Z, wWithcut ceducing the error coefficient. {€&] The 


compensated system is shown in Figure 11. Tc keep the error 
coefficiert ccnstant, the compensator used was a sirple lag 
network (7,< 7, ). 
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FIGURE 11 Lag Compensated System 


A ccupensator which would meet the recuired 
specificaticrs was calculated by conventional methods as a 
check or tke fregram's performance. The Bode piot cf the 
uncompensated and conventionally compensated system is shcwn 
in Figure 12. The values of +7, and +r, fer Figure 11 which 
would meet the design requirements were determined tc he 
ei, = 1C., 7, = 100. The compensated system's respcnse 
using these values for the compensator is shewn in Figure 
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1.2 


1a, 


THAOUT 


FIGURE 13 


ZO 4.0 6.0 8.0 OF 
Time (sec) 


Corventionally Ccmpensated System Respense 


37 


The coupfensated system was redrawn as Shcwn ans 
Figure 14 using the standard program transfer tlocks 
available for system Simulation. 


FIGURE 14 CADS Lag Compensated Systen 


stardard second crder response of 6 = ,3, 
W = 0.75 was chesen as the desired response for the 
n 
Optimizaticn frogram to match. The limits placed cn the 


cptimizaticn frogram were .001 < a <1 and 20 15< ae <1 
Arkitrary vaiues tc begin optimization were specified as 
FE = .017 anc Z = .1 in the logarithmic centers of the 


Co Co 
search 2zcnes. CAES determined 7, = 8.078 and T, = 67.47 


as the optints farameter settings. Figure 15 shows the 
desired respense and the program ccmpensated systen 


response. 


The respense cf the system compensated by CADS 1s more 
nearly the desired response than is the conventicnaliy 
ccupensated syster. One should remember that the specified 
response is for a true second order system whereas, the 


compensated system is forth order. Therefore, a perfect 
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match cf desired versus actual responses could nct he 


cktained. 
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FIGGRE 1& cCADS Compensated Syster Response 
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rc. CASCACE LEAL COMFENSATICN 


The ccmplexity of the next preblem to ke solved Ey CADS 
was extendec to five free parameters. The plant shcwn in 
figure 1€ was to be used tec follcw a unit amplitude 
Sine-wave input of 200 rad/sec. The output amplitude sas to 
re almcst exactly the same as the input amplitude, ané the 
output cculd not lag the input by mcre than 109. Two 


secticns cf cascade compensation were to be used. [8] 


— ae 


—— 


FIGURE 10 Plant for Lead Comrensation 


The specified requirements may be interpreted as an open 


lecp gain > 15 dk énd a phase angle ~ 90° at W = 200 from 
n 
@ Nichols j;lict. A cut and try soluticr on a écde flot 


shewed tkat a gain of 3 XK 106 and ttso phase lead 
comrensatcrs with a double zero at Z = 70 and a double pcle 
at P = 7€C will satisfy the ciosed loop magnitude and fhase 
reguirements. Figure 17 shows the system response and 
desired respcnse obtained for these values. The sagnitude 
of the ccampensated system's respense is 93% of the desired 


response and legs ry 8.12°. 
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The ercklem was then run on the optimization prcgran 
with the kleck connections as shown in Figure 18. The free 
parameters were the poles and zeros of the compensators and 
the gain of the plant. fhe desired data CULVE of 
XDATA = f€1r(z00t) was generated from the standard second 


order step response by setting 46 = 0, WH = 200 and using 
n 


X(z) as tke desired response. The initial searck zone 


Jinits were specified as 600 < P  <800, 60 < Z << 80, 
2 i 


2.€ X 106 < G(3) < 3.2 xX 106. 


FIGURE 18 CACS Block Diagram for Lead Compensaticn 


The optinizaticn program sclution went to the lower 
litits fcr Ecth poles and the upper limits for roth the 
zeros and tke gain. The limits of the search zones were 


relaxed te 500 < Pp. <600, 80 < Z << 9C and 
a at 


3.2 X 106 < €(3) < 3.4 X¥ 106. The program again places the 


free variakties cn the limits of P. = 500, 2 . = 9C, and 
ni i 


G(3) = 3.4 X¥ 108. The system and desired response for these 
values is shcewn in Figure 19. The system magnituce and 
Ehase are much closer to the desired cutput response than 


tke respcnse cf the conventionally designed ccmpensator. 
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Miesphase Gifference is only 4.339. Continued relaxation of 
the kouncaries produced the compensated system shcwn in 
Figure 2C. Figure 21 is a plot of the system's response for 
these values. The system's response is improved in that it 
lags the irput Ey oniy 3.469 and the magnitude 1s 
essentially tre same as the input magnitude. 


sinQoot-4#") 


FIGUEE 20 CALS Compensated System Elock Diagran 


The opftignizaticn program was activated at T = Q 
although tke prcktlem specifications were for the steady 
state resrcnseé. The problem was rerun with the optimization 
precess started after the transient had died out. The same 
values fcr the cptinun compensator were Ccktained. 
Apparently the szajll initial transient did not effect the 
prceklem scluticn. 
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EF. SERVC SYSTEM COFPENSATIGN 


The fcregcing examples of the Simulaticn - optimizaticn 
prcegram were simple examples which could obviously be 
guickly sclved with the standard cut and try metheés. A 
mcre complex and challenging example cf pregram usace is 
resented ky the system Shown in Figure 22. 


The system is an operaticnal servo drive mechanism with 
multiple inputs and discrete level feedback. This bighly 
ncenlinear system's output was to follow the input as closely 
as possikie and in steady state (t2 75ms) there was tec he 
very little reise ripple. The free parameters of the systen 


are the fcles anced zeros of the compensatcr, ow and the 


roles cf the ncise suppressor C . 
n 


Tc sisulete and optimize the system it was redrawn using 
the availabie program blocks as shown in Figure 23. Lupe ng 
Situlaticr it was found that the current limiter for as was 


net needed Lecause Ee < 25 A and the limiter was removed 


tc decreasé frogram run time. The desired response {XDATA) 
waS writter as a set of three eguaticns. These equations 
were then used to replace the second order step respense 
eguaticns in the standard program. DRVIN(1), (4) and _ (5) 
were alsc writter as a set of eguaticns and placed in 
Subroutine ELANT. The discrete level feedback to bicck two 
was achieved by making DRVIN(2) = INTGER(THA(12)). These 
changes tc tte mair program and FLANT were all that were 
neccessary tec simulate this system. The insplementation of 


these changes to the program is shcwn on pages 80 and 81. 
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The system cptimizaticn was initially Eroken intc two 
parts. This was te reduce the number of variables that were 
to be optimized per run and to obtain near cptimum starting 
values fcr the free farameters. The above steps were taken 
in an effort to decrease the computer time required for a 


scluticn. Tke first run was tc cptimize the compensator, 
E¢ independent of the noise supfresscr. The seccnd 


Cptimizaticn cun was to select values for the noise 


SUPELESSCL, Cu The rationale behind this separaticn was 
that the two circuits perform different functicns and should 


therefcre be initially separable in their effects cn the 
System. The final optimization run was to Fe made with all 
free paragreters available to the program fcr optimization. 


Tke search zcne centered on the values found above. 


Figure z4 shows the response for the syste as 


Criginally ccmpenseted. The oftimizaticn trun to set the 


values fcr e resulted in Z = 43.5, Pe = 21.0, 
Z| = 47.5, EF = 5$2.C. Figure 25 shows the simulaticn of 
the system tuSing these values. The initial velocity 


overshoot has keen reduced and the average velocity after 
the trarsient aprears more equally distributed akove and 


kelow the desired velocity. 


50 


ssuodsoeyq maqsks oaras Teutbht3zo hz auNnoTA 


(des) sUTL 
Cts 0 OT*O 80°0O 90°00 . p0°O c0°O °0 


eae ea ee ea ae a a ee A a ee a ee 
a —— SS a Ne ee ea ee Se eye Se ee 


WJ 
© 
none © 


Senet 


0°0¢ 
(DeSs/uTtT) AATOOTSA 


OV 


51 


T 
asuodsay ua3yskg 9 p3szesuadmoy saqyd Sz TNNOTA 


(99S) owt, 


ae Ot 70 9020-9020 70200 Z0°0 


oc 01 


Ss 


Os 
(des/uTt) AAQTOOTAA 


O°Ot 


0° OP 


fhe initial cptimization run for setting the valves of 


tke noise suppresscr used the same cost function whick had 


been used for all previous optimization runs 

J = f ext dt. This resulted in ¢g¢ = 0.2 (lcwer bound), and 
W = 15CO (upper tkound). Figure 26 shcws the result of 
D 


using this cost function was to further decrease the 


cvershoot. Ecwever, it allowed larger switching cr noise 
transients as a resuit of weighting large transient errors 
mcre than the lesser noise jitter. Tc overcome this 
prctlem, tke ccst function was changed to J -f [Err] *t*dt. 
This was tc weight the steady state errcrs mscre heavily than 
tke transient errcrs. Using this cost function the values 


set by tke cptinization prcgram were Sf = .z and W = 1430. 
n 
The damping factcr, d , was again placed on the lower limit. 


The initial cvyershoot was still imprcved cver the original 
System's respcense Eut the switching transients were not 
improved cver the previous optimizaticn trial. The results 
of the sinuletion run using these values is shown in Figure 
7 
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A final optimization run was made with all six 
variables free. The performance index was changed tc the 
weighted cost function shown in equation V. The weichting 
factcr was ccmputed to equally weight the start up transient 
and the steacy state responses in an effort to rceduce the 


Switching transients. 


J =f/{Errjdt 0.0 < T < 0.045 

J = f4.1*[Err{i dt 0.045 <t < TF (V) 
The optinta ccmpensatcr values established by CADS were 
P. = 7A ea Z = 52.5, PS = 52125 Z = 48.1, 6 = 0.22 and 
a = 126€. Figure 28 shows the fully compensated system's 
response. The initial overshoot has been reduced and the 


average velccity in steady state is clcser tc the desired 
value tkan the criginal system's response. However, the 
transient errer cn shut off is still present. Rak Liem ie. 
Summarizes tke results of the optimizaticn of the servo 


systen. 


The CEU time required for the optimizaticn precess was 
net decreased Ly spelitting the problem intc two separate 
parts. The time required fcr each optimization run cn the 
individual ccrrensators was the same as the time reguired tc 


optimize tke ccmpliete system with six variatles. 
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This freblem has presented an excellent example cf the 
necessity tc carefully select the cost function which will 
Geasure the systea's performance. Defining a perfcrmance 
index which will weight the nore objecticnable 
characteristics of a system's response so that they are 
eliginatec or reduced beccmes difficult as the system's 
ccuplexity ircreases. The performance index, J =/(Err)2 dt, 
was not adequate in its treatment of the ncise and switching 
trarsients fcr the above prcblen. The cther perfcrmance 
indexes tsed were also marginal in their effects upon 
Farameter oftisization. Although the ccst functicr was 
reduced to a sathematically ccrrect finipun, the systen's 
perfcrmarnce tas not the best that could ce achieved. The 
system's perfcrmarce waS only optimum due tc the defirition 
of the ccst fturcticn. The system was simulated usinc the 
ccmréensatcr values determined from the last cptimizaticn run 
with the excerpticn of 6 which was increased to Ci Oe 
This simulaticn was made based upon engineering judgment of 
the effects cf varying dé . Figure 29 shows that this change 
in the damping factcr produced a system response which was 


nearer tke desired response than any of the frevious rtns. 
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pa LISCUSSICN 


The okjective of this thesis was to investigate the 
feasibility cf develofing a computer program which weuld 
optimize a variety of control systems' respcnses in the time 
dcmain. ‘he fregram develcped during the investication 
preves that time dcmain optimizaticn cf contrcl system 
responses is net only feasible but desirakle due tc the 
readily interpretable results of the optimization pfrecess. 
CaCS requires appreximately 200K bytes cf ccuputer core when 
the high speed printer plot is used for graphical display of 
the output ard 230K when outputting calcomp gLapas. These 
ccre reguirenents are a maximum and could te reduced Fy more 


careful, prefessional programming. 


The ccuputer time required for CALS to arrive at a 
scluticn is déefendent upon 

(1) the crder of the system being simulated, 

(2) the area cf the search zone determined by the 

upper and lower bounds cn the variable farameters 

and 

(3) the nearness of the starting guess to the 

optinup parameter values. 
Every trial value cf a parameter selected by FOXPLX requires 
a complete sirulation of the system in crder to evaluate the 
system's response and compare it with the desired resfecnse. 
The program run time is therefore, T,., = the number of 
trials xX system simulaticn time. A high crder system may 
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reguire twenty seccnds of CPU time for Sipnulation. If 300 
trials are reguired to determine the cptimum farameter 


values, tke tctal CPU time weuld be 100 minutes. 


Example Lower Starting Upper Number CPU Time 


Problem Bound Guess Bound one Required 
EL XL XS XU Trials 
P -300 34.3 =e 55 7min03sec 
(pseudo) 
-O01 ~O1 0.1 PAPAS. 24min53sec 
Cc Ou api 1.0 225 24min53sec 
400 450 500 
400 450 500 
d 90 95 100 2,000 68min 
90 95 100 


3.4x10© 3.5106  3.5x106 


10 20 : 30 
40 50 60 
e 400 500 600 230 4nr 
3 a 220 
eS O 1300 1400 
Table III 


fable III presents a summary of the search zcnes and 
tines reguired to cktain a solution for scme cf the frecblems 
ccnsidered iz this thesis. The amount of CFU time reguired 
for a soluticn, e¢specially fcr problem JIII-E, may seen 
EXCEessive. Ecwever, there are several consideraticns which 
BHeauld Ee mace pricr tc arriving at this ccnclusion. 

(1) The equaticns of the system do not have to be 


written, programmed nor dektugged if the system's 
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compcenent transfer functions are known. 

(2) Time demain reguirements do not have to be 

translated irto frequency domain specificaticns 

for system sinulation and design. 

(3) A systematic search is carried cut to oktair 

the cptinum farameter settirgs. This assures that 

with valid -krounds on the variables an acceptakle 

scluticn will be obtained with the fips t 

Optimizaticn attempt. 
The time reguired tc perform the above steps in the éesign 
precess Ey ccnventional means may result if wany more hours 
of CPU tine than if the program CADS were used fren the 


beginning of the design process. 


Several means of reducing the computer time recuired 
fcr cptinizaticn were previously outlined in secticn II. 
The most Significant reduction is obtained Ey keepine the 
nugker cf integrations reguired for system sSimulaticn to a 
Mirimum. Elcck diagram reduction of the system shculd be 
acccmplisked whenever possible. The example of the dard 
Lecnard Grive system skhcwn on page 26 can ke reduced tc the 


Sinple system shown in Figure 30 by blceck diagram reduction. 


FIGURE 3G Elock Reduced Ward~Lecnard Drive 
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Section IJI.E. showed the effect of using different 
cost functicns tc measure the quality cf the cptizized 
response. Alttough the integral error squared criteria is 
often used tc judge a system's performance, the user shculd 
carefully ccrsider how to best define a cost functicn which 
will preperly penalize deviations of the system response 
from the desired response. A small error (err < 1) =scuared 
keccmes even Smaller. If all the errors are small, the IES 
is a valid ccst functicn but if the system response invclves 
larce and sg#all errors a weighted cost functicn will have to 
be used. One grethcd cf arriving at a properly weighted cost 
furction is te simulate the system using first estimates of 
tke variakle farameters and recording the sur of the eérrers 
over the different portions of the respense. A ratio cf the 
errers can tken be used tc arrive at proper weichting 


factors fcr each time section cf the respcense. 


BOXELX will ccntinue the optimization frecess, wcrking 
in the seventh significant digit until it can no lcnger 
reduce tke ccst function. Often an acceptarle soluticn for 
the system farameters has been found lorg bkefore the ccst 
function kas teen reduced to its minimum. A judicious use 
cf CADS gay ke made by evaluating the system response after 
ten to fifteer minutes of run time to see if an accertable 


ecluticn has Leen found. 


B. SONCLESICNS 


Time domain cptinization using the CADS program is a 
simple, straight forward process which dces not require an 
in-derpth analysis cf the system being optimized. Eccncmic 
use of CPU times does dictate that intelligent starting 
values and bcunds ke placed on the variable parameters which 


are to be optinizec. 


The frogram 1S a readily usable tool for simulation 
without crtinsization. It is ccmpetitive with, if not 
Superior to, other common simulation routines when 
Situlatirng typical control systems. This feature alcne is 


expected to Iring the program into commcn usage. 


Cc. RECCHMENTATIONS FOR FUTURE WORK 


(1) All integral calculations performed during this 
investigaticn were done in dcuble precisicn for increased 
accuracy cf the system respcense. The possibility cf using 
Single precisicn calculations should be investigated tc 


decrease ccre reguirenents. 


{Z) The akility to begin the optimization process at some 
time greeter than zero shculd be provided. This will 
necessitate nucdificaticn of the data input cards and block 


€guations sc that initial conditions can be entered. 


(2) The frograg presently requires that external fcrcing 
furcticne (CEVINS) and the variables which are to be 
optimized Fe specified by Fortran IV statements flaced 
within tke krEcdy cf the pregram. A method cf reading these 
specificaticrs frcem data card input should be develcped sc 
that the vwuser will not have to "shuffle" cards in the 


prcgram deck. 


(4) The craphical output of CADS was all that was necessary 
for the investigaticns conducted in this thesis [Fut a 
prcevisicn fcr numerical output of selected responses shculd 


ke previded fcr detailed analysis. 


(5) The feasibility of reducing the numker of significant 
digits BCXPLX uses should be studied as a seans of recucing 
cptimizaticn tine. Also some criteria might ke develcrfead to 
re-ect a system's response before TF is reached if it is 


deterfined tc te unacceptable. 
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(6) A method of automatically relaxing the boundaries on the 
variables being optimized when they go to their limits 


should be developed. 


(7) The standard cost function provided and all user 
developed cost functions should be normalized. This would 
permit a more direct and easier comparison of a system's 
"goodness" when several different integration step sizes or 


run times have been used in the optimization process. 
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APPENDIX A 


Block Data Card Format 


11 114 120 440 #460 
ELKCCD=EEVV G P Z 

CC = PCSITICN NUMBER OF THE BLOCK 
[ = TYEE OF BLOCK (NUMBER) 

EE = INECI NCDE NUMBER 

VV = OUTEDT NODE NUMBER 

G = VAIUE CF GAIN 

— = VALUE CF THE POLE* 

Z = WAICE CF THE ZERO* 


* For kLlieck type 4, § is read into the P locaticn and 


Wn is read into the Z pasition. 

* Fer kElock type 5, Sh is read intc the P EOSTTION 
and A is read into the Z position. 

* Fer Elecck tyfe 6, ye is read inte the P position 


and © is read into the Z position. 


EXAHFELES 
G 
PLKC 11=0102 10. 
G 
BLKC22=0504 1. 26 
G ie 
BIKC23=0405 1. ys 24 
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BLK 1iC4=1006 


BLK115=1207 


BLK 126=01708 


*See Table I. 
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APPENDIX B 
Plant Flow Chart 
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FEAC INFLT CATA 
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